Intrinsic Dynamical Fluctuation Assisted Symmetry Breaking in Adiabatic Following 
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Classical adiabatic invariants in actual adiabatic processes possess intrinsic dynamical fluctua- 
tions. The magnitude of such intrinsic fluctuations is often thought to be negligible. This widely 
believed physical picture is contested here. For adiabatic following of a moving stable fixed-point 
solution facing a pitchfork bifurcation, we show that intrinsic dynamical fluctuations in an adiabatic 
process can assist in a deterministic and robust selection between two symmetry-connected fixed- 
point solutions, irrespective of the rate of change of adiabatic parameters. Using a classical model 
Hamiltonian also relevant to a two-mode quantum system, we further demonstrate the formation of 
an adiabatic hysteresis loop in purely Hamiltonian mechanics and the generation of a Berry phase 
via changing one single-valued parameter only. 

PACS numbers: 03.65.Vf, 05.40.-a, 05.45.-a, 37.10.Gh, 45.20.Jj 



Introduction - Adiabatic theorem is about the dynam- 
ical behavior of a Hamiltonian system whose parameters 
are changing slowly with time. It constitutes a funda- 
mental topic in Hamiltonian mechanics [l}. For example, 
Einstein was among the first to recognize the importance 
of classical adiabatic invariants in understanding quan- 
tization 0]. In recent years, there are still considerable 
interests in several aspects of adiabatic theorem in both 
quantum mechanics 0] and classical mechanics [3] • 

Classical adiabatic following is the subject of this 
study, but our findings are also relevant to certain quan- 
tum systems. We start from the fact that classical adia- 
batic theorem is not an exact theorem: an actual adiabat- 
ically evolving trajectory fluctuates around an idealized 
solution predicted by the adiabatic theorem. Adiabatic 
invariants hence possess intrinsic dynamical fluctuations 
(IDF's) The magnitude of such fluctuations, typi- 

cally proportional to the rate of change of adiabatic pa- 
rameters, becomes extremely small in truly slow adia- 
batic processes. So except for special quantities that can 
accumulate IDF during an adiabatic process [4, 8] (e.g., 
in calculations of dynamical angles), IDF does not seem 
to be interesting or physically relevant. As shown in this 
Letter via both theory and computational examples, this 
perception is about to change. 

Bifurcation phenomenon is ubiquitous in nonlinear sys- 
tems and it is of fundamental interest to many top- 
ics in physics, among which we mention localization- 
delocalization phase transitions and symmetry breaking 
[§- 14 1 . Here we consider the adiabatic following of a sta- 
ble fixed-point solution, which, as a result of a varying 



(a) 



R 





Fixed point: 
Where should I go? 



stable fixed point 

unstable fixed point 

FIG. 1: A supercritical pitchfork bifurcation in a classical 
Hamiltonian system parameterized by R. The solid (dashed) 
lines represent the phase space locations of stable (unstable) 
fixed points as a function of R. (a) When R exceeds Ri, 
a stable fixed point splits into two stable and one unstable 
fixed points. One then wonders which branch the system will 
land on as R adiabatically increases beyond Hi. (b) The 
bifurcated fixed points merge back when R exceeds Ri. If R 
first increases from R < Ri to R > R2 and later returns to 
R < Ri, then an adiabatic hysteresis loop may be formed. 



adiabatic parameter, moves towards a supercritical pitch- 
fork bifurcation. As schematically shown in Fig. 1(a), 
two new stable fixed points and one unstable fixed point 
emerge after the adiabatic parameter (denoted R) slowly 
passes the bifurcation point at R = R\ . It is then curious 
to know among the three fixed points, which fixed-point 
solution the system will land on and whether the selec- 
tion is predictable. It is found, both theoretically and 
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numerically, that IDF is crucial for a deterministic and 
robust selection between the symmetry-connected pair of 
stable fixed-point solutions, regardless of how slow the 
adiabatic process is. As such, through crossing the bi- 
furcation point, a tiny IDF is amplified to a macroscopic 
level after the bifurcation: it determines the fate of the 
trajectory afterwards by "forcing" the system to make a 
selection between symmetry-breaking solutions. We term 
this as deterministic symmetry breaking, because there 
is no need for external noise to initiate the symmetry- 
breaking. The selection process is robust because, un- 
like symmetry-breaking selection processes studied pre- 
viously [3], it is independent of the dynamical details. 
Figure 1(b) depicts an interesting situation that involves 
a second bifurcation at R = R2, after which the three 
fixed-point solutions merge back to one stable fixed point. 
We shall show that this case may allows us to generate a 
Berry phase by manipulating one single- valued parame- 
ter only. 

Preliminaries - We recently developed a general de- 
scription of IDF in classically integrable systems 
which can be reduced to a rather simple form for sta- 
ble fixed-point solutions in phase space. In particular, 
let us consider a one-dimensional Hamiltonian H(q,p; R), 
with (q,p) being the canonical variables, R a system pa- 
rameter to be tuned slowly, and its stable fixed point 
solutions denoted by [q(R),p(R)]. According to the tra- 
ditional picture offered by classical adiabatic theorem, a 
stable phase space fixed point has zero action, so when 
R varies slowly, the system must retain its zero action 
as an adiabatic invariant and therefore must follow the 
instantaneous fixed point [q(R),p(R)]. This is however 
a picture without IDF. The actual time evolving state 
(q,p) deviates from [q(R),p(R)}, 
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q = q + Sq, p = p + Sp. 



(1) 



where (6q, Sp) are IDF on top of the idealized adiabatic 
solution [q(R),p(R)] with R — R(t). It is straightfor- 
ward to see why (Sq, Sp) has to be nonzero: were it in- 
deed zero, then by definition of a fixed point, we have 

dH^M = 9H(q,pM) = 0j indicating that the current 

state (q,p) cannot evolve and hence the system can never 
do adiabatic following with a moving fixed-point solu- 
tion. This indicates that nonzero (Sq, Sp) is not due to 
nonadiabaticity. Rather, it is intrinsic and must exist for 
adiabatic following to occur. Our theory applied to this 
case [i| gives (see also [U) 
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FIG. 2: (color online) Impact of intrinsic dynamical fluctu- 
ations on the crossing of a pitchfork bifurcation, (a) The 
system's actual state is on the left of the instantaneous fixed 
point, (b) Immediately after the bifurcation, the actual state 
is still shifted to the left of all the fixed points, (c) As two 
stable fixed points split further, the actual state gets trapped 
by, and starts to adiabatically follow, the stable fixed-point 
solution on the left. 



and (•) denotes an average over all possible initial con- 
ditions of (5q,Sp). As seen from Eq. ([2]), so long as 
exists, then the scale of IDF is proportional to 4? = V, 
the speed of adiabatic manipulation. The magnitude of 
(Sq, Sp) is then deceptively small (for a sufficiently small 
V), but nonzero in general. 

We now examine how nonzero (Sq, Sp) impacts the adi- 
abatic following of [q(R),p(R)} that eventually undergoes 
a supercritical pitchfork bifurcation aX R = R\. First, 
right before the adiabatic parameter R reaches the bi- 
furcation point Ri, the system's actual state (q,p) must 
deviate from the instantaneous solutions [q(R),p(R)] by 
(Sq, Sp). Without loss of generality and based on Eq. (j2J, 
the actual time-evolving state (q,p) is assumed to be 
slightly shifted to the left side of the instantaneous fixed 
point [q(R),p(R)]. This is illustrated in Fig. 2(a), where 
the periodic orbits (associated with a fixed R) around 
[q(R),p(R)} are also shown. The second stage is illus- 
trated by Fig. 2(b), where R only slightly exceeds R\. 
There the bifurcation already occurs but the actual state 
does not "feel" the bifurcation yet: it continues to stay 
on the left side of all the three new fixed points. So the 
actual left-shifted state is located on an orbit (if R were 
fixed) surrounding all the fixed points (this is confirmed 
in our numerical studies). The line shown in Fig. 2(b) 
passing through the unstable fixed point represents the 
separatrix. In the last stage, R further increases, the 
two stable fixed points split further, with the stable fixed 
point moving to the left capturing the actual state [see 
Fig. 2(c)]. During the ensuing adiabatic following, the 
actual state then adiabatically follows the instantaneous 
fixed point on the left. Clearly then, if, as R approaches 
a bifurcation point, IDF can induce a definite shift (to 
the left or the right) of the actual state with respect to 
[q(R),p(R)\, then the system will be trapped, determin- 
istically, by one of the two stable fixed-point solutions 
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after the bifurcation. Note that, exactly because of IDF, 
the actual state always stays away from the vicinity of 
the unstable fixed point (the extremely slow part of a 
separatrix). As a result the diverging time scale associ- 
ated with the whole separatrix does not affect adiabatic 
following here. This understanding will be confirmed in 
our following numerical experiments. 

Model Hamiltonian - Here we turn to a concrete Hamil- 
tonian system with supercritical pitchfork bifurcations. 
Specifically, we choose 



H = --q 2 - Ry/l - q 2 cos(p) + Ay/l - q 2 sm(p) (4) 

in dimensionless units, with c > A > 0. Inter- 
estingly, this system has two bifurcation points at 
R x = -Vc 2 -A 2 and R 2 = Vc 2 - A 2 . Define 77 = 
y/l - (R 2 + A 2 )/c 2 and fi = arctan(-A/i?). In the 
negative regime of R, the stable fixed point is (q,p) = 
(0, fi — 7r) for R < R\, which then bifurcates into two sta- 
ble fixed points at (^fr), fx — 7r). In the positive regime of 
R, the two stable fixed points are at (^77, //) for R < R2 
and then merge back to one stable fixed point at (0, p). 
Note that this model Hamiltonian is invariant under the 
joint operation of space reflection (V) q —> —q and time 
reversal (T) t — > —t. Therefore, for cases with only one 
stable fixed point, the solution itself has the same VT 
symmetry, and for cases with two stable fixed points, the 
pair are transformed to each other under the VT opera- 
tions, with each individual fixed point being a symmetry- 
breaking solution. The phase space locations of these 
stable fixed-point solutions for various fixed values of R 
are shown by the solid lines in Fig. 3. 

Another motivation to choose H in Eq. (|4]) is that it 
describes a two-mode many-body quantum system on 
the mean-field level. Consider the following mean-field 
Hamiltonian in dimensionless units (K = 1), 



-M 2 ) 



-R-iA 



-R + iA -c(\b\ 2 - \a\ : 



(5) 



where a and b are quantum amplitudes on two modes, c 
represents the self-interaction strength, and R ± iA rep- 
resents inter-mode coupling. By rewriting a = \a\e l ^ a , 
b = \b\e"t' b , p = 4>b - 4> a , and q = \b\ 2 - \a\ 2 , the time 
evolution of q and p, as obtained from the Schrddinger 
equation for this quantum model, becomes precisely that 
under the classical model Hamiltonian H in Eq. (j4j . The 
fixed points of H become eigenstates of H m ; the adiabatic 
following when crossing a bifurcation for H is mapped to 
the issue of adiabatic following as degenerate eigenstates 
of H m emerge. One possible realization of H m is a Bose- 
Einstcin condensate (BEC) in a double-well potential, 
with the imaginary coupling constant A implemented via 
phase imprinting on one well 15|. 7f m may also be re- 
alized in nonlinear optics by using two nonlinear optical 
waveguides with biharmonic longitudinal modulation of 
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FIG. 3: (color online) Solid lines depict phase space loca- 
tions of fixed-point solutions for our model Hamiltonian in 
Eq. (|4)| as values of R are scanned. The solid and empty dots 
represent the actual states of the system in numerical ex- 
periments of an adiabatic process, where R is first increased 
from 7? = —0.25 < Ri to R = 0.25 > R2, and then returns 
to R — —0.25, at a constant speed V = 10 -6 . Long ar- 
rows indicate the moving direction of the actual states as R 
changes. The four double-well plots schematically show the 
symmetry-breaking consequence in the context of a two-mode 
BEC Hamiltonian. c = 0.2, A = 0.1. Variables plotted here 
and in Fig. 4 are in dimensionless units. 



the refractive index [l6(. So our detailed results below 
are relevant to both classical and quantum physics. 

Theory and numerical experiments - Applying the the- 
ory of IDF to the Hamiltonian in Eq. one obtains 



(Sp) = and (Sq) 



dR 



for 



VR 2 +A' 2 R 2 +A 2 -cVR 2 +A 2 dt 

R < Ri [17|. Therefore, as R increases from R < Ri, 
^ = V > 0, and (Sq) is definitely negative. Return- 
ing to the phase space plot in Fig. 3, this means that 
if R increases from R < Ri, the actual state (on av- 
erage, to be more precise) is slightly left-shifted from 
the fixed point at (0,/x — 7r). So after passing the bi- 
furcation at R — Ri , the system is expected to adiabati- 
cally follow the left symmetry-breaking solution (q, p) = 
(— 77, /j — 7r) for i?i < R < and (q,p) = (— 77,/j,) for 
< R < i?2- Similar theoretical results show that, if 
R decreases from R > i? 2 , then the actual state must 
be slightly right-shifted from the fixed point at (0,/j), a 
prediction consistent with the above-mentioned VT sym- 
metry of the system. So after crossing the bifurcation 
point at R = R2, the system should adiabatically follow 
the other symmetry-breaking solution (q,p) — (77, /x) for 
< R < R 2 and (q,p) = (77, /x - tt) for Ri < R < 0. 
As shown in Fig. 3 (solid or empty dots), our numer- 
ical results based only on Hamilton's equation of mo- 
tion confirm this prediction. For the results shown we 
have set V = I0~ 6 to ensure a slow process. Indeed, at 
all times the difference between the actual states (dots) 
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and the instantaneous fixed points (solid lines) is invis- 
ible to our naked eyes. Yet, the small IDF does assist 
in a symmetry-breaking choice regarding which of two 
stable fixed-point solutions is adiabatically followed by 
the system. In the language of H m , after the system 
passes the bifurcation at R = Ri owing to an increasing 
R, q = \b\ 2 — \a\ 2 becomes appreciably negative, so one 
mode develops more population than the other, signal- 
ing a clear delocalization-localization transition induced 
by IDF. The opposite population imbalance occurs when 
the system passes the bifurcation at R — R2 with a de- 
creasing R. Furthermore, joining these two manipulation 
steps together so that R returns to its very original value 
in the end (which is already the case in Fig. 3), we clearly 
observe the formation of an adiabatic "hysteresis" loop in 
phase space. That is, by increasing and then decreasing 
R, a navigation loop in phase space is formed because 
the adiabatic following in the forward step and backward 
step lands on different symmetry-breaking branches fl8j j . 

There is one subtle point to be clarified: Our theory of 
IDF [see (Eq. H}] is about quantities (Sq) averaged over 
all initial conditions of (Sq, Sp), but what determines the 
adiabatic following is the actual Sq in a single process. 
The Supplementary Material contains a detailed analysis 
on this point. In particular, we show that if initially 
there is a deviation of Sq from (Sq), then the difference 
Sq — (Sq) oscillates with time and later, as R approaches 
the bifurcation point, this deviation becomes negligible 
as compared with (Sq). It is for this reason that the 
definite sign of (Sq) is equivalent to the definite sign of 
Sq, which hence justifies our theory based on (Sq). 

Berry phase generation via one single-valued parame- 
ter - As a final interesting concept, we discuss the gener- 
ation of a Berry phase using one single-valued adiabatic 
parameter R. This is made possible by IDF and bifurca- 
tions. In particular, the navigation loop in Fig. 3 shows 
that the time-evolving states of H m trace out a nontriv- 
ial geometry after increasing R from R < R\ to R > R2 
and then returning R to its initial value. Let \ip\eit) and 
fright) be the eigenstates of H m (with i?i < R < R 2 ) 
mapped from the left and right fixed-point solutions of 
H. Assuming exact adiabatic following with the instan- 
taneous adiabatic eigenstates, the Berry phase generated 
along the navigation loop is analytically found to be 
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(6) 



In our numerical experiments, we choose to integrate the 
Berry connection using the actual states during the phys- 
ical process (this simple method will not account for the 
nonlinear Berry phase correction studied in Ref. 0). It is 
found that for the shown regime in Fig. 4, the agreement 
between theory and simulations is excellent, with tiny 
but visible differences. Such visible differences remind us 
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FIG. 4: Analytical and numerical Berry phase results vs c, 
for the quantum system described in Eq. ([5]) with A = 0.1. 
The theoretical results are from Eq. while the numerical 
results are obtained directly by integrating the Berry connec- 
tion along the actual path of adiabatic evolution. 



that a direct numerical integration of the Berry connec- 
tion along the actual time evolution path is not necessar- 
ily reliable because it could accumulate the effect of IDF 
[4j. Nevertheless, the fair agreement shown in Fig. 4 con- 
firms our analytical result, demonstrates the feasibility of 
Berry phase generation using only one single-valued adi- 
abatic parameter, and verifies from another angle that 
IDF is important for understanding adiabatic following 
in the presence of bifurcation. 

Conclusion - Bifurcation greatly amplifies subtle in- 
trinsic fluctuations that are beyond classical adiabatic 
invariants. This leads to a selection rule regarding which 
of two symmetry-connected stable fixed-point solutions 
may be adiabatically followed. In cases of multiple bifur- 
cation points, adiabatic hysteresis loops in phase space 
and the generation of Berry phase by manipulating one 
single-valued parameter are also shown to be possible. 
Our findings are of fundamental interest to both classical 
systems and quantum many-body systems on the mean- 
field level. The implications of this work for symmetry 
breaking in fully quantum many-body systems should be 
a fascinating topic in our future studies. 



Appendix 

In this Appendix we discuss the difference between (Sq) 
and Sq, using the same notation as in the main text. It 
is important to clarify this point because our theory of 
intrinsic dynamical fluctuations (IDF) is about the quan- 
tity (Sq) averaged over all initial conditions of (Sq,Sp), 
but, as indicated in our main text, what determines the 
symmetry-breaking adiabatic following is the actual Sq 
in a single process without the averaging. Specifically, 
we need to show that, before the system crosses the bi- 
furcation point, the (positive or negative) sign of Sq in 
an actual process is always the same as the sign of (Sq) 
predicted by our theory. Without loss of generality we 
choose to discuss the case R < R\ as an example. 
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For our model Hamiltonian, the T matrix evaluated at 
the fixed point is found to have zero diagonal elements. 
Hence we can write it in the following form 



-A 
B 



(7) 



Next we expand the Hamilton's equation of motion to 
the first order of Sq and Sp, we have 



-A 
B 

-A 
B 



Sp 
Sq 

p-p 

g-Q 



(8) 



where the instantaneous fixed point (q, p) is located at 
[0, arctan(— A/R) — ir]. Note that in general the off- 
diagonal elements A and B are i?-dependent. Never- 
theless, to gain physical insights and to develop a simple 
analytical result from the above equation, we consider a 
small time segment during which the T matrix can be 
regarded as a constant, and p moves at a rate M due to 
an increasing R. We then obtain 



dp . 

Tt = - Aq; 

| = Bp-BMt, 



(9) 



a first-order differential equation that can be integrated 
directly. Upon direct integration we find the analytic 
solution 



Dcos(VABt) 



M 
1 



q 

p = -ADsin(VABt) + Alt, (10) 
where D is one integration constant. This then indicates 



M 



5q = Dcos(VABt) - — ; 



Sp 



-ADsm(VABt). 



(11) 



Equation (fTT|) clearly shows that Sq and Sp are oscil- 
lating solutions. In particular, the oscillation amplitude 
in Sq is seen to be the initial difference between Sq and 
the averaged quantity (Sq) — —M/A (for the concerned 
time segment). Even more importantly, as the system 
approaches the bifurcation point at R = the term 
—M/A is the dominating term because \A\ sharply de- 
creases for R approaching R±. This makes it clear that 
the sign of Sq must agree perfectly with the sign of (Sq). 

To verify this theoretical understanding, we performed 
numerical experiments and quantitatively compare in 
Fig. [5] the numerically found Sq with (Sq) obtained from 
the theory of IDF (see the main text), for an increasing 
R (R < Ri). As seen in Fig. [SJa), the actual numer- 
ical Sq indeed oscillates around theoretical (Sq). As R 
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FIG. 5: (color online) Comparison between numerical results 
of Sq and theoretical (Sq) obtained from theory of IDF. Sys- 
tem parameters are c = 0.2, A = 0.1, and ^ = 10~ 6 . (a) 
Initially the canonical variables (q,p) are set exactly on the 
instantaneous fixed point. The initial Sq is hence zero, which 
differs from the theoretical nonzero (Sq). (b) Initially Sp and 
Sq are precisely chosen according to the theory of IDF (see 
the main text). In this case, Sq is seen to agree with (Sq) at 
all times before the bifurcation, which is about to occur at 
R = —0.173 for parameters used here. 



approaches closer to i?i = —0.173 (beyond the shown 
regime of R in Fig. [5]), the oscillations of Sq around (Sq) 
remain small and yet (Sq) becomes more and more neg- 
ative. Therefore, the sign of Sq is always identical with 
the sign of (Sq) when the system is close to bifurcation 
crossing. As a second confirmation of our insights above, 
Fig. EJb) shows that, if initially we set Sq = (Sq) and 
Sp = (Sp), then the ensuing time dependence of Sq as 
found from our numerical calculation agrees exactly with 
that from the theory of IDF. 

We finally note that, the ultimate reason why our an- 
alytical solution above for (Sq, Sp) during a small time 
segment is already so useful is again linked with adia- 
batic following. That is, since R changes slowly, the off- 
diagonal elements of T and M will all change slowly, and 
hence the equation of motion for Sq, or for the deviation 
Sq — (Sq) , is expected to adiabatically follow the solution 
given in Eq. (fTTj) . This is interesting because it suggests 
a theory of higher-order fluctuations in the intrinsic fluc- 
tuations. 
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